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Abstract 

Human mobility is a key factor in spatial disease dynamics and related phenomena. In computational models 
host mobility is typically modeled by diffusion in space or on metapolulation networks. Alternatively an effective 
force of infection across distance has been introduced to capture spatial dispersal implicitly. Both approaches do 
not account for important aspects of natural human mobility diffusion does not capture the high degree of pre- 
dictability in natural human mobility patters, e.g. the high percentage of return movements to individuals' base 
location, the effective force of infection approach assumes immediate equilibrium with respect to dispersal. These 
conditions are typically not met in natural scenarios. We investigate an epidemiological model that explicitly cap- 
tures natural individual mobility patterns. We systematically investigate generic dynamical features of the model 
on regular lattices as well as metapopulation networks and show that generally the model exhibits significant dy- 
namical differences in comparison to ordinary diffusion and effective force of infection models. For instance, the 
natural human mobility model exhibits a saturation of wave front speeds and a novel type of invasion threshold 
that is a function of the return rate in mobility patterns. In the light of these new findings and with the availability 
of precise and pervasive data on human mobility our approach provides a framework for a more sophisticated 
modeling of spatial disease dynamics. 
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I. INTRODUCTION 

The 2009 outbreak of a novel subtype (H1N1) of influenza A and its subsequent worldwide spread, the 
recent emergence of new human infectious diseases such as SARS in 2003, and the recurrent seasonal 
ourbreaks of influenza epidemics illustrate the growing importance of understanding the dynamics of 
human infectious diseases (IH3). Key to understanding spatial dynamics in particular is an accurate as- 
sessments of human mobility patterns as infectious diseases spread among different locations due to 
movements of their host. Despite recent advances EJ |5) comprehensive data on mobility is typically 
unavailable, and modelers have to make reasonable assumptions when implementing host mobility 
in models. Often it is assumed that hosts move randomly (Fig|T][a)) in the system yielding reaction- 
diffusion dynamics EJ ISlfTTl . An alternative heuristic approach captures spatial dynamics without ex- 
plicitely accounting for host dispersal [12j [13) • Instead an effective force of infection between spatially 
separated populations mimics the effect of disease transmission across distance. Typically, this force is 
assumed to be proportional to the prevalence of the disease in one of the locations. Yet, because this 
approach lacks the explicit incorporation of host movements, a systematic analysis of conditions under 
which it is applicable is difficult. 




Figure 1: Models for human mobility: Patches and arrows represent individual populations and travel flux, re- 
spectively, (a) Diffusive dispersal: indistinguishable individuals travel randomly between different locations gov- 
erned by the set of transition rates o) nm . (b) Natural human mobility: Individuals labeled k travel from their base 
location k to connected locations m and back with travel rates o)^ k and respectively. 

Recently, human movement patterns became accessible based on pervasive mobility proxies EJ |5j 
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fTiHITl . One of the key findings of these studies confirmed the intuitive notion that humans spend most 
of their time in small sets of particular locations (home, work, shopping sites, etc.) and a person's mo- 
bility occurs predominantly between these individual-specific locations. Furthermore, a typical charac- 
teristic is the existence of one or two locations that function as an individuals base location, e.g. their 
homes to which individuals typically return before they travel to another place. A key feature of human 
mobility therefore is a bi-directional pattern in their trajectories among small sets of salient locations 
contrasting diffusion processes in which individual agents eventually visit every location in the entire 
system. Topologically, natural human mobility patterns can be described by individual mobility net- 
works that possess a hub-and-spokes structure, in which a central hub represents a base location and 
a limited set of places connected by spokes the set of popular destination locations. Spreading phe- 
nomena across a large spatial scale occurs by virtue of interactions of agents that possess overlapping 
individual mobility networks. 

In preceding papers and in this conference we have presented a stochastic model that explicitly ac- 
counts for the natural human mobility patterns described above [14j[T8]. In particular, it respects the 
fact that individuals typically return to their unique base location before they travel to a new destina- 
tion (recently a similar approach has been used by Balcan and Vespignani 02)). Here we investigate 
properties of this model and demonstrate that the dynamic consequences of natural mobility patterns 
are profound. Its basic aspects are depicted in Fig. [I] In this model mobility of the entire population is 
represented by a set of overlapping individual mobility networks. In the language of complex network 
theory each individual mobility pattern consists of a central node (the base location) connected to a set 
of accessible destinations (connected nodes) in the aforementioned hub-spokes topology. 

Although mathematical metapopulation models have been proposed that are able to capture natural 
human mobility patterns (20j [21] it has remained elusive to what extent and under which conditions 
such models exhibit dynamic features that are qualitatively different from ordinary reaction- diffusion 
processes. It is unclear how these systems can be related to paradigmatic reaction- diffusion systems. 

In the present article we show that the dynamics exhibits profound differences, indeed, as compared 
to ordinary reaction- diffusion systems. We concentrate on the analysis of epidemics on regular lat- 
tice and complex metapopulation networks. On lattices we obtain a generalization of the paradigmatic 
Fisher-Kolmogorov equation that describes wave propagation in reaction-diffusion systems. We show 
that the spatially continuous version of our model exhibits travelling wave solutions and compute their 
front velocities as a function of system parameters. Contrary to reaction- diffusion systems that exhibit a 
monotonic and unbounded increase of the front velocity with increasing travel rate, our model predicts 
an upper bound for front velocities. We show that the front shape strikingly differs from those predicted 
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by reaction-diffusion systems and is more robust in response to changes in parameters. We introduce a 
commuting ratio parameter, a quantity present only in the natural mobility model, and investigate the 
front velocity as a function of it. Analysing a fully stochastic system in regular lattices as well as complex 
metapopulation networks we find that a global outbreak of a disease is determined by a novel threshold 
that is determined by the typical time spent away from individuals' base locations. 

II. NATURAL HUMAN MOBILITY AND DISEASE DYNAMICS 
A. Disease Dynamics on a Metapopulation 

We consider a system of populations labeled m = 1,...,M and assume that in each population an 
epidemic outbreak can be described by a compartmental SIR-model, i.e. 

Im + S m ► 2I m , I m ► R m , (1) 

in which the reactions govern infection due to the interaction of infected (J) with susceptible (S) individ- 
uals at rate a, and recovery of an infected individuals at rate /3, respectively. The number of individuals 
in a population is given by N m = S m + I m + R m . The spread of an epidemic across the set of M populations 
is governed by the exchange of individuals between populations. The most prominent and conceptually 
clearest ansatz is diffusive dispersal between populations in which individuals of each class move from 
location n to m at rate co mn , i.e 

X n ► X m , (2) 

where X n represents I n > S n or R n . The rates o) mn generate an equilibrium distribution N n of individuals 
among M populations. Assuming that for a pair of populations exchange rates are nonzero, detailed 
balance is fullfilled, i.e. N n /N m = co nm /a) mn . In the following we will assume that the system is in equi- 
librium with respect to dispersal, i.e. N m = N m , yielding the following mean-field dynamical equations: 

d t I n = aS n I n /N n - pi n + {(D nm I m -a) mn I n ) 
dtSn — ~ocS n I n IN n + ^ {(O nm S m — OJ mn S n ) 
Rn — N n — S n — I n . (3) 

Note that in the reaction- diffusion system individuals are indistinguishable apart from their infection 
status and move about randomly between all available locations m. This approach has been employed 
both in complex networks of coupled populations as well as simplified lattice models fTT1l22l . 
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The relation to spatially continuous models is best illustrated in a system of linearly aligned popula- 
tions separated by distance /, at locations x = nl, uniform population size N n = N = const, and travel 
between neighboring populations only, i.e. co nm = cod n -i >m + to8 n+ i )m . The overall uniform rate to is re- 
lated to the waiting time in a given location t = co~ l . In the limit /, t — ► with D = I 2 It this model yields 
the 1-d reaction- diffusion system 

d t j = ajs-/3j + Dd 2 x j 

d t s = -ajs + Dd 2 s, (4) 

where j{x,i)l = I n IN n and s{x,t)l = S n /N n . These equations are related to the Fisher-Kolmogorov 
equation [8j|9). For sufficiently localized initial conditions j{x, t = 0) this system exhibits travelling waves 
with front velocity 



c = 2l^(a-/3)/T~Va), (5) 

that monotonically increases with the global mobility rate to. 

In order to account for individual mobility networks that exhibit base locations and natural recurrent 
movements we propose the following generalization of Eqs. (T),(2): We assume that individuals can be 
grouped into subpopulations defined by two indices, n and k. The first index determines the current 
location n> the second the base location k. Generally, the dispersal dynamics is then governed by a set 
of reactions: 

X n V X w (6) 



0) 



This implies that individuals of class k possess their specific dispersal rate matrix co^n that is condi- 
tioned on the base location k. The rate ix)\ in determines how individuals of type k travel from location n 
to m, for fixed k the matrix co^n represents the aforementioned individual mobility network for individ- 
uals of type k. The dynamical system, incorporating disease dynamics, is given by 



where l\ and S\ are the number of infecteds and of susceptibles of type k located at n> respectively. N n 
denotes the total number of individuals in location n> i.e. N n = LfcC^ + + ^n)- Note that if the rates 
(o„ m are independent of k } we recover the ordinary reaction- diffusion case. 

In the following we consider the case of overlapping hub-spokes networks corresponding to com- 
muting between base and destination locations. This imposes restrictions on the rates (i)^ mt it implies 
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that co k m = if k ^ n and k ■£ m. That means individuals of type k that are located at m cannot travel 
to n without returning to k first. We further assume that o) k mk = o)~ , i.e. the return rate is uniform for all 
k and m. This assumption implies that individuals typically spent the same amount of time in distant 
locations before returning to their base. Assuming that travel of the entire system is equilibrated, we 
obtain 

k 1 + Lm^ k mk lu k km 

where N n is the stationary number of individuals located in population n and N k = Y. n N k is the total 
number of individuals of type k (i.e. they belong to base location k). 

An important limiting case is a situation in which mobility rates are large compared to the rates as- 
sociated with the infection and recovery dynamics, i.e. co^^co - » a, /3. In this case detailed balance is 
fulfilled for infecteds and susceptibles separately and the last terms in Eq. (7) vanish. If we assume that 
o) k mk l(x)~ « 1 which implies that individuals belonging to k remain at their base most of time, Eq. (7) can 
be reduced to the effective force of infection model ITS! : 

^-I k = ccS k j:e km I m -pi k , (8) 

where I k = X m is the number of infected individuals belonging to location k and coupling strengths 

TYl 171 

£nk - Zm Pn^P^IN are explicitly related to travel rates and p™ = N n IN n is the occupation probability. 
Hence direct coupling represents a special case of our model, see also EU . 

In order to investigate the dynamic consequences of natural human mobility patterns as captured by 
Eqs. |7| ,we consider a system analogous to the one-dimensional spatially homogeneous system lead- 
ing to the reaction-diffusion Eq. f4). We consider a 1-d lattice of populations of size N separated by a 
distance /, assume next-neighbor coupling and allow only infecteds to travel (relaxing this restriction 
does not change the main results but eases the analysis). We denote the number of infecteds at their 
base location by I" and the number of infecteds at neighboring locations [n - 1) and [n + 1) by J~ or J+, 
respectively. This yields 

In p In O) 

where S n denotes the number of susceptibles at n. o) + and o)~ denote forward and return rates, re- 
spectively. In the corresponding dynamical system we can approximate S n ,I^ by their continuous 
counterparts: J* +1 — ► ^{x ± I) « J* + IV I* + f-AZ*. In an equilibrated homogeneous lattice the size 
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Figure 2: Front propagation in natural human mobility models. Left: Front velocity c(o)) as a function of travel 
rate co for a reaction -diffusion and the natural mobility model. Numerical results of the stochastic simulations with 
N = 10 4 agents per site are depicted by symbols for reaction- diffusion (blue) and natural mobility (red) models. 
Analytical results for the reaction- diffusion system, Eqs. (S) and natural mobility model, Eqs. (TT) are depicted by 
dash- dotted and dashed lines, respectively. Right: Front velocity c{e) as a function of the commuter ratio e (for 
0) = 1). The dashed curve illustrates the analytical result, Eq. (15} . Crosses and circles depict results of numerical 
solutions of the mean-field equations (7) and stochastic simulations of the natural mobility in a spatial SI model, 
respectively. 

of a single population remains constant during an epidemic. Introducing concentrations u n = I n IN, 
Vn = (/+ + I~)/N and w n = (J+ - I~)/N this yields: 

d t u - a(\-u-v){u+v + D/\v + l\/w) + (x)~v-2(x)*u 

d t v - 2oj + u-oj~v (10) 
d t w = -a)~ w, 

where D = I 2 12. The third equation is solved by w - 1 and for t » l/co~ yields w « 0. We can 
therefore discard w leaving only the first two equations in (10). Steady states are u = 0, v = and u = 
(x)~ /{2o) + + a)~), v = 2o) + /{2o) + + (x)~). In the non-zero steady state concentration of infecteds in one city 
is given by u+ v = 1. In order to compare this system to the reaction- diffusion model, we calibrate 
both systems such as to keep the total flux of individuals between two particular locations equal in both 
systems. 
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III. RESULTS 



In epidemiology key questions concern conditions under which an epidemic can spread. In this case 
it is generally the first task to compute or estimate the speed at which an epidemic proliferates through- 
out the entire system. On regular lattices with dispersal only among adjacent site this task is equivalent 
to computing asymptotic wavefront speeds. In complex network topologies alternative quantities are 
useful, e.g. the time to reach the epidemic peak. Below, we investigate the velocity of epidemic propa- 
gation on a lattice, and subsequently provide results in a complex network topology. 



A. Front velocity in lattice systems 

Using the traveling wave ansatz /(x, t) = g(x - ct) for u and v> performing a linear stability analysis 
of the disease-free state, we find that the system defined by Eqs. (10) exhibits traveling wave solutions 
with front velocity given by 



c= ——. (11) 

a + oj~ +2o) + 

If the forward and return rates co + and co~ are significantly different, two extreme cases can be consid- 
ered. In the limit w + ^0we find c —> as expected, i.e. no propagation can be sustained in the limit of 
individuals not leaving their base. If the backward rate co~ is small, the system is determined exclusively 
by the forward rate co + . 

It is instructive to first consider a balanced system, i.e. co + = co~ = co. In this case Eq. ([TTJ simplifies to 

2\/6Da(0 

c= — . (12) 

a + 3o) 

The front velocity as a function of infection rate a and travel rate to as well as results of stochastic nu- 
merical simulations are depicted in Fig. [2] For comparison, the front velocity of the reaction- diffusion 
scenario with the same global travel rate co is depicted as well. The velocity in this case is given by Eq. J5) 
which increases with travel rate according to - y/co. In contrast, the natural mobility model exhibits a 
saturation of the front velocity with increasing travel rate. From Eq. (12) it follows that the asymptotic 



value of the velocity is proportional to the reaction rate lim^oo c = 2aV2D/3. 

The existence of a saturation is a consequence of natural mobility patterns that are restricted to in- 
dividual mobility networks. Likewise the unbounded increase in front velocity in a reaction- diffusion 
systems is a consequence of the unnatural assumption that increasing the travel rate also increases an 
individual's access to the entire system. In the more realistic natural mobility model increasing co only 
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increases the rate of travel between the base and the two neighboring sites and does not imply faster 
coverage of the entire system. 

In contrast to the reaction-diffusion system with only one rate parameter to, the natural human mo- 
bility model possesses two travel rates, co + and co~ . The total flux between two neighboring locations is 
given by F nm = (o + N% + co~ N^, for m = n ± 1. In equilibrium, detailed flux balance requires F nm = F mn . 
Comparing to the reaction-diffusion system, the total flux is given by F® m = coN. In order to compare the 
dynamics of both systems quantitatively it is plausible to gauge both systems such that flux is identical, 
i.e. 



Q) + N% + Q)-N% l = Q)N. (13) 

To simplify the analysis it is convenient to introduce a commuting ratio £ = co + /co~. In a situation in 
which individuals dwell at their base location most of the time, the commuting ratio is a small, £«1. 
We can express forward and return rates in terms of the global travel rate co and the commuting ratio £ 
according to: 

,1 + 26: , l + 2£ 

co = co and co = co. (14) 

2 2£ 



With these definitions expression (TT) can be rewritten as 

2cc£V2D (2+1/e) 

c = -. (15) 

l + 2£ + 2a£(l + 2£)- l o)- 1 

Fig. [3] depicts the front velocity as a function of the commuting ratio at the fixed travel rate {to = 1) and 
illustrates a significant change of c as a function of £. Since the commuting ratio is not defined for 
ordinary reaction- diffusion systems, this effect cannot be captured in these systems. 

A surprising result in the behavior of the front velocity for very low travel rates is depicted in Fig. [2] 
which illustrates a noticeable deviation of the stochastic system from both, the analytical prediction of 
Eq. Q~T) and the results of the numerical solution to Eqs. (7). The deviation of the numerical solution 
from the Monte Carlo simulations for small co is due to the finite number of agents per site. We observe 
a crossover from a linear scaling with co (symbols) towards the numerical mean-field solution (solid 
blue line). The regime of low travel rates effectively corresponds to high infection rates {a » co). This 
implies that an outbreak takes place almost instantaneously in a neighboring location and and epidemic 
essentially jumps from one location to the next. The rate of hopping is proportional to Nco, i.e. to the 
flux of individuals between locations, where N is a typical of individuals per site. The crossover from 
discrete to continuous behavior occurs when a - coN, i.e. co c - aIN. Note that the slow convergence 
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Figure 3: Front speed in the limit of small travel rates. Bottom Left: Front velocity c(oj) in the regime of low 
travel rates co. Symbols reflect stochastic simulations for different values of particles per site. Dashed and dash- 
dotted lines represent analytical results given by Eqs. (TT) and (5), respectively Solid blue line corresponds to the 
numerical solution of Eqs. (7) and the solid black line represents the analytical result of Eq. (19} . 
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Figure 4: Front shape at the leading edge: Front shape for different global travel rates oj in the reaction- diffusion 
system (crosses) and natural mobility model (circles) combined with local SIR dynamics. 



of the velocity towards zero with decreasing travel rate can be understood qualitatively: Let's consider 
just two locations with agents that can travel between them. Without loss of generality we consider an 
SI epidemic. At the beginning of the epidemic the number of infecteds in the second location is small 
and we can linearize the standart SI dynamics for the second population, i.e. 

dj 2 



dt 



«a/2 + <U7i, 



(16) 



where we neglect the backward flux of the individuals from the second location. Integrating Eq.[l6|by 
means of the integrating factor jzit) = e at J Q r dr j\ (r) and using the solution of the SI model for the first 
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Figure 5: A novel type of invasion threshold. Left: The panel depicts the results of stochastic simulations of an SIR 
epidemic in combination with natural mobility patterns. The attack ratio p(o)~) as a function of the return travel 
rate co~ exhibits a critical transition and vanishes for return rates exceeding a critical value. The function p{o)~) is 
shown for a 1-d lattices system with 100 sites with N agents/ site. Insets: (a) p{o)~) for a scale-free network (y = 1.5) 
of 10 3 nodes populated uniformly with (TV) = 250, k m [ n = 5 and /c max = 50; (b) p(o)~) for an Erdos-Renyi network 
with 500 nodes and (/c) = 10. Results were averaged over 50 realizations. Global travel rate was kept constant at 
0) = 1. Epidemic parameters are a = 1, /5 = 0.1. Right: The attack ration as a function of both travel rate parameters 
w ± for a lattice system. Solid lines represent constant global rate co curves with logarithmic spacing. The total 
flux increases from bottom left to top right. Green regions corresponds to sustained outbreaks, red reflects the 
extinction regime. Epidemic rates, lattice and averaging parameters are the same as in the left panel. 



location y'i(r) = 1/(1 + ae~ at ) with a = 1 - 7i(0)/y'i(0) yields 

l + ae~ at 
1 + a 



j 2 {t) = e at a) In 



(17) 



As t — ► oo we have ln(l + ae~ at ) ~ 0, and thus 72(f) ~ (oe at . The time lag between outbreaks in both 
populations is given by At(^) = Tziq) - T\{q) y where T\ and T2 are times when the concentrations of 
infecteds attains some threshold q in the first and second location respectively. It follows that 



q q 

Ar(g) - In In In a, 

aa) 1 - q 



(18) 



and for q = 1/2 we obtain for the velocity c - At 1 and thus 

c~ ("In— -lnaj . (19) 

This expression reproduces the slow convergence of the front velocity towards zero with decreasing 
travel rate which is in agreement with the results obtained by numerical solution of the mean-field prob- 
lem. 
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Figure 6: Effects in stochastic system: Left: The behavior of the attack ratio p as a function of the ratio of return 
rate and oj~ global rate oj. Other parameters are identical to those of the 1-d system in Fig. (5). Right: Time lag 
t{o)~) as a function of co~ for the natural mobility model (circles) and ordinary reaction- diffusion (dashed line) on 
a scale-free network (with scaling exponent y = 1.5 and 10 3 nodes with an average of (N) = 250 individuals per 
site. The curves depiced averages over 50 stochastic realizations. 

B. Front shape 

Front velocity is tightly connected with the shape of the front at the leading edge. In Fig. [2] front 
shapes as obtained from numerical solutions for both reaction-diffusion and natural mobility SIR mod- 
els are shown. We observe that the reaction-diffusion slope is much more strongly affected by the global 
travel rate co. This is in contrast to the system with natural mobility patterns in which, just like the front 
velocity, the front shape converges to a fixed state as co increases. The particular front shape is related to 
Kendall waves that have been empirically observed |23] . 

C. Invasion Thresholds in natural human mobility models 

In addition to questions concerning the velocity of disease propagation, it is of fundamental impor- 
tance to assess the conditions under which an epidemic propagates at all. Usually such a condition takes 
the form of a threshold in a system parameter. Thea most prominent example is the basic reproduction 
number, given by Ro = a//3 for a SIR model. It quantifies the a number of secondary cases caused by a 
single infected individual in a totally susceptible population [TJ. If R > 1 an outbreak occurs, otherwise 
the epidemic wanes. Another threshold parameter in the metapopulation reaction- diffusion framework 
is the global invasion threshold. It represents the minimal required flux of individuals traveling between 
two locations Q~T] in order for a disease to propagate spatially. 
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One of the most striking properties of the natural mobility model of Eqs. {7J is the existence of a novel 
type of threshold that is only determined by the return rate 0)~ , or equivalently by the typical time an 
individual spends at a distant location. The existence of this threshold is evident from Fig. [5] that depicts 
the attack ratio p (the total fraction of infecteds during an epidemic) as a function of the return rate co~ 
on a) a one-dimensional lattice, b) and Erdos-Renyi network and c) an uncorrelated scale-free network. 
For low return rates the attack ratio is close to unity, as expected, a global outbreak occurs. However, 
with growing values of the return rate, the attack ratio drops almost to zero, i.e. no global outbreak 
occurs. The regime of high return rates corresponds to small dwelling times at distant locations. This 
implies that an infected does not have sufficient time to transfer the disease to susceptibles in unaffected 
locations before returning home. This effect is absent in reaction- diffusion systems. This novel type of 
threshold is a direct consequence of the properties of natural human mobility. 

To assess the mutual impact of all travel parameters, i.e. the forward and return travel rates as well 
as the total flux on the dynamics we calculated the attack ratio for various parameter values on a homo- 
geneous lattice. The results are presented in Fig. [5] Note that our model exhibits the global ^-limited 
invasion threshold that is also present in ordinary reaction-diffusion systems. However, for large total 
flux co, the system exhibits a global outbreak only if the return rate co~ is sufficiently small. Increasing 
the return rate the system enters a region that lacks a global outbreak. Consequently, only the return 
rates o)~ is a limiting factor. This novel threshold can be estimated analytically. In the same spirit as 
introduced recently 1 18) [22] we find the following threshold relation: 

No) 

(Ub - 1) > 1. (20) 

p + ti)~ 

The inverse of the sum of return rate o)~ and recovery rate /3 provides the typical time an infected in- 
dividual spends on a distant location in the infected status. Using the relation co = 2o) + o)~ / {2o) + + o)~) 
(compare |14)), we can write explicitly 

0)~ ti) + 

— <—2N{Ro-l)-l. (21) 

Note, that from (20) the empirically observed scaling p = p(Na)/to~) follows, see Fig. [5] and [6j The 
figures exhibit the expected collapse of the data according to this scaling. The pronounced difference 
between ordinary reaction- diffusion systems and the natural mobility model is also captured in Fig. [6] 
that illustrates the epidemic peak time r - f I{t)tdtl f I{t)dt on the return rate. The figure shows that 
by varying the return rate one can substantially increase the peak time in the natural mobility model, in 
contrast to the reaction-diffusion model. 
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Figure 7: Natural human mobility and disease dynamics in realistic scenarios. (A): The minimal spanning tree 
of a multiscale mobility network in the United States. This network represents the skeleton of the most important 
transportation routes in a more complete network 0]|24) • (B) , (C) : Comparison of the evolution of the SIR epidemic 
with initial outbreak in Los Angeles for the natural mobility model (B) and ordinary reaction- diffusion model (C). 
Snapshot of the number of infecteds after approximately 4 weeks. Color encode from high (red) to low (blue) the 
relative number of infecteds per county. 

IV. DISCUSSION 

An unprecedented amount of information on human mobility available today requires adequate 
models that correctly capture the key features of natural human mobility in order to understand, de- 
scribe and predict the dynamics of human mediated contagion phenomena. In the present work we 
pursued this goal and formulated an approach that can account for important features of natural hu- 
man mobility that are intuitive and observed empirically. The model is based on a metapopulation 
approach assuming well-mixed local populations and explicit incorporation of individual mobility pat- 
terns conditional on base or home location. We considered regular bi-directional movements of the host 
between base locations and accessible destinations, systematically analyzed the model and compared 
it to established modeling approaches, e.g. effective force of infection and reaction- diffusion systems. 
We demonstrated that the latter are limiting cases of the natural human mobility model at low and high 
travel rates, respectively. For a regular lattice we derived a generalization of the Fisher-Kolmogorov equa- 
tion and found that contrary to the reaction-diffusion approach, the front velocity of the epidemic does 
not increase unboundedly with increasing global travel rate, but saturates at a maximum level. 

Although results for lattice and artificial random network topologies are extremely helpful in gaining 
fundamental insight into the dynamics and consequences of natural human mobility patterns on the 
patterns of disease spread, they at best mimic real world scenarios. A future task will be to investigate 
to what level these novel effects prevail in more realistic settings, i.e. real world mobility networks on 
which individuals move and transmit disease. We hypothesize that effects that are so dominant in par- 
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simonious systems of the type described by the model of Eqs.[7]will also be present in more complex 
settings. Evidence for this has recently been revealed in a multiscale metapopulation system [19). In 
order to illustrate the pronounced difference in disease dynamic patterns that are generated by ordinary 
reaction diffusion models on one hand and natural human mobility models on the other, we simulated 
both systems on the backbone of a realistic, multiscale mobility network in a real geographic setting. 
Fig. [7] illustrates snapshots of the timecourse of disease spread generated by both models. Nodes in the 
network are approx. 3000 counties in continental United States and links between them resemble the 
traffic flux. We observe a significantly smaller spreading speed in the natural mobility scenario com- 
pared with the reaction- diffusion model. This implies that estimates of spreading speeds, also in these 
more realistic setting, could have been overestimated in the past by models that rely on ordinary diffu- 
sion as a dispersal mechanism. 

Although the study of natural human mobility on disease dynamics and related human mediated 
contagion processes requires more attention in future investigations, the results presented here as well 
as in our previous work [14j [18) will serve as a useful guide for developing more reliable large scale 
computational models for disease dynamics. The substantial differences of ordinary reaction- diffusion 
dynamics and the novel type of natural human mobility disperal suggest that disperal mechanism are 
among the most important modeling ingredients that require particular care when implemented in large 
scale computational models that are designed to make quantitative forecasts. 

D.B. andV.B. acknowledge support from the Volkswagen Foundation. D.B. acknowledge support from 
EU-FP7 grant Epiwork. 
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